arXiv:hep-ph/021131 lvl 20 Nov 2002 


UWThPh-2002-34 


Two-pion contribution 
to the muon magnetic moment 1 

G. Ecker 

Institut fur Theoretische Physik, Universitat Wien 
Boltzmanngasse 5, A-1090 Wien, Austria 

The two-pion contribution to hadronic vacuum polarization can be extracted from 
r decay data when isospin violating and radiative corrections are taken into ac¬ 
count. When the dominant corrections are applied to the photon-inclusive decay 
r~ —* !/77r - 7r°[7], one obtains a shift A a p = (—12.0 ±2.6) x 1CP 10 for the anoma¬ 
lous magnetic moment of the muon. The shift appears to be too small to reconcile 
the determinations of hadronic vacuum polarization from existing t and e + e~ 
data. The reliability of electromagnetic corrections in the photon-inclusive r de¬ 
cay is examined. 


1 Status of g fl — 2 

The magnetic moment of a particle with spin S p , 

€-Tl —* 

9p = 9p 2^~c Sp ’ (1) 

is expressed in terms of its gyromagnetic factor g p . With polarized muons in a prop¬ 
erly tuned storage ring [1], the spin precession frequency is directly proportional to 
the anomalous magnetic moment a p = (g p — 2)/2. On the basis of the most recent 
experimental result from Brookhaven [2], the present world average is 

a® xp = (11659203 ± 8) x lO -10 . (2) 


The standard model prediction for consists of three parts: 

„SM _ QED _|_ weak , had 

a /j ~ a /j, ' 

The first two contributions are known with very high accuracy [3]: 

a QEE> = (H658470.6 ± 0.3) x 10“ 10 

a w eak = ( 15.2 ± 0 . 1 ) x 10" 10 . 
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Fig. 1: Contribution of lowest-order hadronic vacuum polarization to a 


To set the stage, we may compare the experimental result with the combined electroweak 
contribution, 


a axp - a^ ED+weak = (717 ± 8) x 10~ 10 


( 5 ) 


a difference of 90 standard deviations. 

Although the primary motivation of the Brookhaven experiment was to test the 
weak contribution a™ eak , the main issue at present for a meaningful comparison between 
theory and experiment is to understand the hadronic contribution a kad : 


had had,LO . had,HO _i_ ~had,LBL 

a H — a fi ' a fi ' a Li 


( 6 ) 



8±4 


— 10 . 0 ± 0.6 


By far the most important hadronic contribution is due to vacuum polarization at 
lowest order in a (shown in Fig. 1): 


had,LO _ hvp 


( 7 ) 


The most recent value for the higher-order hadronic contribution [4] is displayed in 


Eq. (6) (all values for are given in units of 10 10 from now on). The hadronic light- 
by-light contribution a h ad - LBL s ^j]] carr j es a l ar ge theoretical uncertainty but at least 


the sign is now established [5, 6]. 

Most of the hadronic vacuum polarization contribution to originates at rather 
low energies (about 70 % for t < 0.8 GeV 2 ). Therefore, nonperturbative methods are 
needed together with experimental input. The most recent evaluation, using either only 
e + e“ data or including r decay data (for the two- and four-pion channels), finds [7] 


hvp _ 


684.7 ± 7.0 
701.9 ±6.2 



( 8 ) 


r 



Anomalous magnetic moment 


3 


The two determinations of aj) vp in Eq. (8) differ by more than four standard deviations 
[7]. This discrepancy could be of experimental origin or the theoretical analysis might 
be inaccurate or incomplete. It is the main purpose of this talk to discuss whether the 
discrepancy could be due to an underestimate of isospin violation. 


2 Hadronic vacuum polarization and isospin viola¬ 
tion 

The contribution of hadronic vacuum polarization at 0(a 2 ) to a M (Fig. 1) is given by 

[ 8 ] 

/»oo 

aj) vp = / dtK(t)oo(e + e~ —> hadrons)(t) (10) 

JiMl 

with a smooth kernel K{t) concentrated at low energies. I discuss here only the two-pion 
contribution that accounts for 73 % of a]) vp . 

A 1 

The two-pion contribution can also be extracted from r decay. In the isospin limit, 

, , _ , _ w . , . . dT(r” —> 7r°7r ~v T ) . „. 

cr 0 (e + e -> 7T + 7T )(t) = h(t) --- (11) 

at 

with a known kinematic function h[t). At the level of accuracy needed to match the 
present experimental precision, a systematic account of isospin violation including elec¬ 
tromagnetic corrections is required. I report here on a recent analysis of those correc¬ 
tions [9, 10]. 

We expect the size of isospin violating corrections to lie somewhere between 


M 2 


= 2 x 10 


-3 


and 


M 2 , — M 2 c 

7T+ _ 7T_ 

MS 


= 0.067 


In a first step one integrates out all heavy fields with masses > m T . This generates 
an electroweak short-distance correction factor for semihadronic r decays [11] Sew — 
1.0194 (in the MS scheme). The second and final step is then to calculate the isospin 
violating corrections in the theory with light fields only. 

The CVC relation (11) gets modified in the presence of isospin violation: 


cr 0 (t) = h(t) 


dT(r —> 7r°7r v T ) Risit) 
dt. Sew 


with an isospin breaking correction function 


Risit) 


1 

G EM (t) 


F v{t) 

f+(t) 


( 12 ) 


2 


(13) 
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The phase space correction factor [12] 




1 + 


3(M*+ - Mgg) 

i - 4Af2 


+ 0[(M 2 + 


Mi 


(14) 


is especially important near threshold. The ratio of form factors (Fy(t) in e + e~ anni¬ 
hilation, f+(t) in r decay) is mainly sensitive to p — to mixing and (to a lesser extent) 
to the width difference P p + — T p o [7, 10]. 

The main task is the calculation of radiative corrections that go into the function 
Gem(£)- As usual, the radiative corrections consist of two parts and only the sum is 
infrared finite and well defined: the exclusive rate with one-loop corrections [9] and 
the radiative rate. I concentrate here on the calculation of the radiative rate [10], i.e., 
r(r“ —► n°TT~iy r "f) under ALEPH conditions [13], with photons of all energies included. 

To describe this decay, we have used a gauge invariant chiral resonance model with 
the following features [10]: 

• Low’s theorem (leading and subleading terms) is manifestly satisfied in terms of 
an explicit representation for the pion form factor f+(t) [14]. 

• The amplitude exhibits the correct low-energy behaviour to 0(p 4 ). 

• The low-energy amplitude is extended into the resonance region using the standard 
chiral resonance Lagrangian [15]. The implicit assumption is that p and (to a 
lesser extent) a\ exchange, which contribute already at 0(p 4 ), are the dominant 
mechanisms at all accessible energies. 

For photon energies < 100 MeV (in the r rest frame), the rate is dominated by 
bremsstralrlung (leading Low approximation). However, under ALEPH conditions with 
all photons included, the bremsstrahlung approximation is not sufficient. The infrared 
finite sum of loop-corrected and radiative rate translates into the function G'em^) shown 
in Fig. 2. 

We are now ready to calculate the total shift in a p vp due to isospin violation in the 
two-pion channel: 


A a^ vp 

H- 


/4M2 


dtK ( t ) 


h(t)- 


dT „ 


r H 


dt 


x (^> -i 
\ -JEW 


(15) 


In Table 1 the various contributions and the total shifts are displayed for two values 
of t ma x- electroweak short-distance correction (5 ew)> threshold correction (KIN), ra¬ 
diative corrections (EM) and form factor ratio (FF). Obviously, the low-energy region 
below 1 GeV 2 is dominating. 

I have only listed our error estimate for the form factor ratio because the uncer¬ 
tainty associated with some electromagnetic low-energy constants appearing in Gem(^) 
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Fig. 2: Electromagnetic correction function Gem(£ ): the full curve corresponds to the com¬ 
plete radiative amplitude, the dashed curve is based on the leading Low approximation 
(bremsstrahlung) for the radiative mode. 


tmax (GeV 2 ) 

Sew 

KIN 

EM 

FF 

total 

1 

- 9.5 

- 7.5 

- 1.1 

6.1 ± 2.6 

- 11.9 ± 2.6 

3 

- 9.7 

- 7.5 

- 1.0 

6.1 ± 2.6 

- 12.0 ± 2.6 


Table 1 Contributions to Aa^ vp from various sources of isospin violation (in units of 10 10 ) 
for two different values of t max (in units of GeV 2 ; f max < m 2 ). 


is much smaller [10]. The short-distance and kinematic corrections are model indepen¬ 
dent (although there could be higher-order corrections in Sew)- The form factor ratio 
is dominated by p — co mixing taken directly from the most recent experimental analysis 
[! 6 ]- 

The total isospin violating correction goes in the right direction towards reconciling 
the t with the e + e~ data but the shift seems to be too small in absolute magnitude. 
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While the signs and magnitudes of the shifts caused by the short-distance correction, 
the kinematic threshold effect and the form factor ratio are well understood, the small 
electromagnetic shift deserves further discussion. 

A first observation is that loop and radiative contributions tend to interfere destruc¬ 
tively (for reasonable infrared cutoffs because only the sum is infrared finite). A more 
instructive exercise is the comparison with the radiative corrections for the inclusive 
rate T(t _ — » duv vpy]) at the quark level, calculated some time ago by Braaten and Li 
[17]. Translating their result into a shift in a^, one finds 

Aa EM,quark = L g , ( 16 ) 

of opposite sign but similar magnitude as in Table 1. There is no fundamental reason 
why the inclusive result at the quark level should agree with the exclusive two-pion 
result. Nevertheless, performing the radiative corrections for the two-pion mode in 
the leading Low approximation for the radiative mode (independent of any resonance 
contributions to the decay amplitude), one obtains accidentally exactly the same value: 

Aa EM,Low = 16 ( 17 ) 

Note that the shift in the inclusive case and in the bremsstrahlung approximation for 
the exclusive channel goes in the “wrong” direction increasing the discrepancy with the 
e + e“ result. 

What is then the origin of the “correct” sign of the shift Aa® M in Table 1? It 
turns out that the difference to the bremsstrahlung value is completely due to the 
subleading terms in the Low expansion for the radiative amplitude proportional to the 
derivative df+(t)/dt of the pion form factor [10]. In other words, the full curve in Fig. 2 
could not be distinguished from the curve based only on the first two terms in the 
Low expansion. To the extent that the pion form factor is known experimentally, the 
shift Aa® M is therefore model independent and certainly independent of details of the 
resonance exchange model for the radiative decay t~ —> t^tx~v t ^. The result depends 
only on the (shape of the) pion form factor. 


3 Conclusions 


The total shift of aj] vp due to isospin violation and radiative corrections in the two-pion 
channel, 

AaJ] vp =-12.0 ±2.6 , (18) 

agrees well with a similar more data-oriented analysis of Davier et al. [7]. Performing 
the shift for aj) vp extracted from the two-pion decay of the r, one arrives at the results 
[7] already displayed in Eq. (9): 


a 


exp 


-«r = { 


33.9 ±11.2 (3.0 a) [e+e~} 

16.7 ±10.7 (1.6 cr) [r] 


(19) 
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Before drawing any far-reaching conclusions about possible evidence for new physics, 
the reason for the discrepancy between the two determinations of aj] vp must be under¬ 
stood. One possible uncertainty has to do with the experimental procedure of applying 
radiative corrections to the e + e _ data where the corrections are much bigger than in 
the t decay. If both the e + e _ result and the raw r data were correct isospin violation 
would have to be more than twice as big as calculated [10]. In view of the discussion 
presented here, especially on the natural size of electromagnetic corrections in the r 
decay to two pions, I consider such a drastic underestimate very unlikely. 

For the time being and pending clarification of the discrepancy between e + e _ and 
r-based extractions of a[) vp , the standard model prediction for the anomalous magnetic 
moment of the muon is in good shape. 


Acknowledgments 

I wish to congratulate G. Martinska, J. Urban, S. Vokal and their team for the excellent 
organization of Hadron Structure ’02. I also thank V. Cirigliano and H. Neufelcl for a 
very pleasant collaboration. 


References 

[1] Y.K. Semertzidis: hep-ph/0211038 and references therein. 

[2] G.W. Bennett et al. (Muon g — 2 Coll.): Phys. Rev. Lett. 89 (2002) 101804. 

[3] E. de Rafael: hep-ph/0208251. 

[4] B. Krause: Phys. Lett. B 390 (1997) 392. 

[5] M. Knecht and A. Nyffeler: Phys. Rev. D 65 (2002) 0730334. 

[6] M. Knecht, A. Nyffeler, M. Perrottet and E. de Rafael: Phys. Rev. Lett. 88 (2002) 071802. 

[7] M. Davier, S. Eidelman, A. Hocker and Z. Zhang: hep-ph/0208177. 

[8] M. Gourdin and E. de Rafael: Nucl. Phys. B 10 (1969) 667. 

[9] V. Cirigliano, G. Ecker and H. Neufeld: Phys. Lett. B 513 (2001) 361. 

[10] V. Cirigliano, G. Ecker and H. Neufeld: JHEP 08 (2002) 002. 

[11] W.J. Marciano and A. Sirlin: Phys. Rev. Lett. 61 (1988) 1815; ibid. 71 (1993) 3629. 

[12] H. Czyz and J.H. Kuhn: Eur. Phys. J. C18 (2001) 497. 

[13] R. Barate et al. (ALEPH): Z. Physik C 76 (1997) 15. 

[14] F. Guerrero and A. Pich: Phys. Lett. B 412 (1997) 382. 

[15] G. Ecker, J. Gasser, A. Pich and E. de Rafael: Nucl. Phys. B 321 (1989) 311. 

[16] R.R. Akhmetshin et al. (CMD-2): Phys. Lett. B 527 (2002) 161. 

[17] E. Braaten and C.S. Li: Phys. Rev. D 42 (1990) 3888. 



